High overtones of Dirac perturbations of a Schwarzschild black hole 
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Abstract 

Using the Frobenius method, we find high overtones of the Dirac quasinormal spectrum for the Schwarzschild 
black hole. At high overtones, the spacing for imaginary part of lo„ is equidistant and equals to Im(_j n +i) — Im(w„) = 
i/8M, (M is the black hole mass), which is twice less than that for fields of integer spin. At high overtones, the 
real part of u> n goes to zero. This supports the suggestion that the expected correspondence between quasinormal 
modes and Barbero-Immirzi parameter in Loop Quantum Gravity is just a numerical coincidence. 
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The Quasinormal Mode (QNM) spectrum is an impor- 
tant characteristic of a black hole. It dominates the late 
time response of a black hole to an external perturbation, 
and, at the same time, does not depend on the way of their 
excitation. Thus, being dependent on black hole parame- 
ters only, the QNMs provide us with the "fingerprints" of 
a black hole, feasible to be seen in the detection of gravi- 
tational waves (for reviews, see [1]). 

The importance of the QN spectrum is not limited 
by the above observational aspects of gravitational waves 
phenomena, since QNMs have interpretation in conformal 
field theory through AdS/CFT and dS/CFT correspon- 
dence [2]. There is also suggestion that the asymptotic 
QNMs of black holes are connected with the so-called 
Barbero-Immirzi parameter in Loop Quantum Gravity, 
which must be fixed to reproduce the Bekenstein-Hawking 
entropy formula within this theory [3] . All this stimulated 
considerable interest to study the QNMs of black holes in 
flat, dS and AdS backgrounds [4]. In particular, the QNMs 
of the Dirac field for different black holes were considered 
in the papers [5], [6], [7], [8]. Nevertheless, the study of 
the Dirac QNMs of four-dimensional black holes were lim- 
ited by low overtones. The low overtones were found for a 
Schwarzschild black hole in [6] , with the help of the third- 
order WKB method [9]. There it was observed that for 
low overtones the real part of the frequency decreases as 
the damping grows. Soon, the analysis was extended to 
the case of Schwarzschild-de Sitter black hole [7] using the 
sixth-order WKB method [10], and also to the case of the 
charged Dirac field in the background of a charged black 
hole [8]. 

In the present paper, we are trying to cover this gap in 
the study of the Dirac QNMs and shall investigate the high 
overtone behavior of the Dirac perturbations of a Schwar- 
zschild black hole. We observed that at high overtones the 
imaginary part of the QN spectrum is equidistant with the 
spacing Im(w„ + i) — Im(w„) = i/8M. It is also shown that 
when the overtone number n is increasing, the real part of 
u) approaches zero. 

The Dirac equation in an arbitrary curved background 
has the form [11]: 

(W^ + r*.) +"*)* = «. (1) 

where m is the mass of the Dirac field, e a M are tetrads, 
are spin connections [11]. The Schwarzschild metric has 
the form: 

dr 2 

ds 2 = f(r)dt 2 - - r 2 (d0 2 + sin 2 9d(f> 2 ), (2) 
f{r) 

where /(r) = 1 — (2M/r), M is the black hole mass. In this 
background, the equation for the massless Dirac field can 
be reduced, after some algebra, to the wave-like equation: 

f^5+^-^(o)*(r*) = 0, (3) 



with the effective potential [12] 

which vanishes at the both boundaries: V(r* = ±oo) = 
0. The parameter /i corresponds to a multipolc number. 
Under the choice of the positive sign for the real part of u> 
(u = ujB,e — iuim, wj{ c > 0), QNMs satisfy the following 
boundary conditions 

*(r*) ~ C± exp(±zwr*), r* — > ±00, (5) 

corresponding to purely in-going waves at the black hole 
event horizon and purely out-going waves at infinity. Then, 
following Leaver [13], we can choose 

*(r*) = exp(iwr*)u(r*), (6) 

where u(r) has a regular singularity at the event horizon 
and is finite at r* — ► 00. 

The appropriate Frobenius series is 

00 

u(r) = f(r) 2s J2 a nf(r) n/2 , (7) 

f(r) s oc cxp(-iwr*), s = -2Miuj, r* — > -00. (8) 

Substituting (6) and (7) in (3) we obtain the five-term 
recurrence relation, which, after a series of Gaussian elim- 
inations can be reduced to the three-term recurrence rela- 
tion 

a n+1 a n (uj) + a n f3 n (u>) + 

Q>n — 1 In 

When u is the quasinormal frequency, the ratio of the 
series coefficients is finite and can be found from the stan- 
dard continued fractions [13]. 

The QNMs are the roots of the inverted continued frac- 
tion: 

> Jn - X 0n-2-a„-37»-2/- /3„+2-a„ + 27~+3/... 

that can be solved numerically as soon as a n (oj), /3„(w), 
7„(w) are found. 

Note, that for the above effective potential, we are not 
able to use the Nollert expansion [14]. That is why in spite 
of the slow convergence of the continued fractions in the 
unmodified Leaver procedure we had to be limited by it. 
This certainly requires much greater computing time to 
perform the computations. 

Low overtones. The low overtones can be obtained 
either with the help of the Frobenius method or, when 
the overtone number n is less than the multipolc number 
/z, applying the WKB formula [9], [10]. Let us remind 
that the potentials with opposite chirality produce the 
same spectrum (see for instance [7] and references therein). 
That is why we shall consider only positive values of u. 
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stratc the following asymptotic behavior: 
Re(w) n « as n — > oo, 



(9) 



Table 1: First ten quasinormal frequencies of Dirac per- 
turbations for /i = 1,2 found by the Frobenius method. 
For modes with n < /i the 6th order (**) and 3th order 
(*) WKB values are given. 

The WKB approach has been used recently when study- 
ing QNMs of different black holes (see [15] and references 
therein) and shown good agreement with numerical com- 
putations for lower overtones. As can be seen from the 
tables above we observe very good agreement between the 
WKB and Frobenius data, and, the higher the WKB or- 
der, the closer the obtained WKB values to accurate nu- 
merical results. 

High overtones. The main difference from what we 
know on the high damping regime for perturbations of 
fields of integer spin (scalar, gravitational, and electro- 
magnetic) is that now the spacing in imaginary part is 
not i/AM, as it takes place for integer spin perturbations, 
but i/8M. The real part of ui (See Fig. 2 and Fig. 3) falls 
down quickly to tiny values which already cannot be found 
with reasonable accuracy by the Frobenius method. This 
occurs already at about n = 400 — 500 for fj, = 1 multipole 
and at about n = 8000 for \i = 2. Note that the greater 
the ratio Imu) n /Reu) n , the more slowly the Frobenius se- 
ries converge, so we had to be careful when finding high 
overtones and to check the convergence by increasing of 
the "length" of the continued fraction until the result for 
uj n will not change. This required, for instance, the length 
of the continued fraction up to 5 millions for n = 7000. 
That is because we were not able to use the Nollert proce- 
dure for the considered effective potential, and therefore, 
a lot of computer time was required to get the convergence 
of the Leaver procedure. That is the reason why we could 
not extend our computation to higher n than that shown 
on Fig 2 and Fig 3. Note also, that for n larger then 500 
it is very time consumingly to compute the modes one by 
one, so we had to "skip" through the hundreds of overtones 
as it is shown on Fig 3. 

Thus, as can be seen from Fig. 2 and 3, QNMs demon- 



Im(w, i+ i) - Im(w„) w -i/8M as n — > oo. (10) 

The formula for the spacing of the imagianary part 
can also be reproduced following [16]. Thus, for highly 
damping modes we can use the Born approximation, where 
the scattering amplitude is given by the formula [16]: 



/+oo 
-OO 



V[r(r*)]e 2lkr dr* 



(11) 



The above integral gets significant contribution only near 
the event horizon. Using the Eq. 4 for V[r(r*)], we find 

S(k) ~ combinations of T(AikM) and T((l/2) + AikM). 

The poles of the amplitude occur when (1/2) + AiMk = 
—n, or AiMk = — n (n > and is integer), i.e. k n = 
in/8M is the required high frequency asymptotic for imag- 
inary part. The same result can be obtained cither using 
the Taylor expansion of the effective potential near the 
event horizon, or without such expansion, but in the lat- 
ter case the expression for the scattering amplitude will 
have a cumbersome form. The singularities of the scatter- 
ing amplitude are shown on Fig 1. 

Conclusion. In this paper we have studied the high 
overtones of the Dirac quasinormal spectrum for a Schwarz- 
schild black hole. The spacing in imaginary part is two 
times less than that for integer spin fields when n — > oo. 
As was shown both numerically and analytically in [17] 
the real part of QNMs for scalar and gravitational per- 
turbations asymptotically approaches a constant equal to 
log3/87r. On the contrary, for the electromagnetic per- 
turbations the real part goes to zero [18]. The real part 
for Dirac perturbations goes to zero, as it happens for 
electromagnetic perturbations. This supports the sugges- 
tion that the expected corresondence between quasinormal 
modes and the Barbero-Immirzi parameter in Loop Quan- 
tum Gravity is just a numerical coincidence [19]. 

The above analysis can easily be extended to the case 
of massive Dirac field [6], [20]. At low overtones, massive 
Dirac pcrtubations [6], similar to massive scalar pertuba- 
tions [21], lead to greater oscillation frequencies and slower 
damping. At high overtones, similarly to the treatment in 
[22] for massive scalar field, one can easily anticipate that 
the massive term will not affect the asymptotic quasinor- 
mal behavior. 
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Figure 1: The absolute value of the scattering amplitude S(k) as a function of k up to a constant factor. M = 1, 



[2] G.T.Horowitz and V.Hubeny, Phys . Rev. D62 
024027 (2000); D. Birmingham, I. Sachs, and S.N. 
Solodukhin, Phys. Rev. Lett. 88, 151301 (2002); 
A.Starinets, Phys.Rev.D66: 124013, 2002; V. Car- 
doso, R. Konoplya and J. P. S. Lemos, Phys. Rev. 
D 68 044024 (2003); E. Abdalla, B. Wang, A. Lima- 
Santos, W.G. Qiu, Phys.Lctt.B538:435-441,2002; 
R.A. Konoplya, Phys.Rev.D66:044009,2002; 
E. Abdalla, K.H.C. Castcllo-Branco, A. Lima-Santos, 
Phys.Rev.D66:104018,2002. 

[3] S.Hod Phys. Rev. Lett. 81 4293 (1998); O. Dreyer, 
gr-qc/0211076. 

[4] E. Bcrti, K.D. Kokkotas, Phys.Rev.D67:064020,2003 
[gr-qc/0301052]; C.Molina, D.Guigno, E. Abdalla 
and A.Saa Phys. Rev. D69 104013 (2004); 
R.A.Konoplya and A.Zhidenko, JHEP (06) 037 
(2004); R.A.Konoplya, Phys. Rev. D 66 084007 
(2002); Phys.Rev.D70:047503,2004; G. Siopsis, 
hep-th/0409262; 

[5] V.Cardoso and J.P.S.Lcmos, Phys. Rev. D63 124015 
(2001); 

[6] H. T. Cho, Phys. Rev. D 68, 024003, 2003. 

[7] A.Zhidenko, Class . Quant . Grav . 21 273 (2004) 

[8] Wei Zhou, Jian-Yang Zhu, Int.J.Mod.Phys.D13:1105- 
1118,2004. 

[9] S. Iyer and C. M. Will, Phys. Rev. D35 3621 (1987) 

[10] R.A.Konoplya, Phys. Rev. D68 024018 (2003) 

[11] N.D. Birrcll and P.C.W. Davies, Quantum Fields in 
Curved Space, Cambridge Univ. Press, Cambridge 
(1982) 

[12] D. R. Brill, J. A. Wheeler, Rev.Mod.Phys. 29:465- 
479,1957; S. Chandrasekhar, The Mathematical The- 
ory of Black Holes (Clarendon, Oxford, 1983). 

[13] E.M. Leaver, Proc.Roy.Soc.Lond.A402:285-298,1985 

[14] Nollert H-P 1993 Phys. Rev. D 47 5253 



[15] R.A. Konoplya, Phys.Rev.D68: 124017,2003; 

[16] T.R.Choudhury and T.Padmanabhan, Phys. Rev. 
D69, 064033 (2004); A.J.M. Medved, D. Martin and 
Matt Visser, Class . Quant . Grav. 21 1393 (2004) 

[17] L.Motl, Adv.Theor.Math.Phys. 6: 1135, 2003: 
E.Berti and K.Kokkotas, Phys.Rev.D68:044027, 2003 

[18] L.Motl and A.Neitzke, Adv. Thcor. Math. Phys. 
7, 307 (2003); V.Cardoso, J.Lcmos and S.Yoshida, 
Phys.Rev.D69:044004, 2004. 

[19] E.Beri, V.Cardoso and S.Yoshida, Phys.Rev. D69: 
124018, 2004: J.Natario and R.Schiappa, hep- 
th/0411267. 

[20] H. T. Cho, and Y-C. Lin gr-qc/0411090 

[21] R.A. Konoplya, Phys.Lett.B550:117-120,2002. 

[22] R.A. Konoplya, A.V. Zhidcnko, gr-qc/0411059. 



4 



Re (w) 



. 175 



0.15 - 



. 125 - 



0.1- \ 



0.075 - 



0.05 



. 025 - 



10 



30 



40 



50 



-Im (w) 



Figure 2: Real part of 10 as a function of imaginary part for fi = 1 multi-pole. 
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Figure 3: Real part of lo as a function of imaginary part for /j = 2 multi-pole. 



